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ABSTRACT 

Objectives This study assesses whether the observed 
declines in HIV prevalence since the beginning of the 
'Avahan' India HIV/AIDS prevention initiative are 
consistent with self-reported increases in condom use by 
female sex workers (FSWs) in two districts of southern 
India, and provides estimates of the fraction of new 
infections averted among FSWs and clients due to 
increases in condom use in commercial sex after 2004. 
Methods A deterministic compartmental model of 
HIV/sexually transmitted infection (STI) transmission 
incorporating heterogeneous sexual behaviour was 
developed, parameterised and fitted using data from two 
districts in Karnataka, India. Three hypotheses of condom 
use among FSWs were tested: (H 0 ), that condom use 
increased in line with reported FSW survey data prior to 
the Avahan initiative but remained constant afterwards; 
(Hi) that condom use increased following the Avahan 
initiative, in accordance with survey data; (H 2 ) that 
condom use increased according to estimates derived 
from condom distribution data. The proportion of fits to 
HIV/STI prevalence data was examined to determine 
which hypothesis was most consistent. 
Results For Mysore 0/36/82.7 fits were identified per 
million parameter sets explored under hypothesis HrjAV 
H 2 , respectively, while for Belgaum 9.7/8.3/0 fits were 
identified. The HIV epidemics in Belgaum and Mysore are 
both declining. In Mysore, increases in condom use 
during commercial sex between 2004 and 2009 may 
have averted 31.2% to 47.4% of new HIV infections in 
FSWs, while in Belgaum it may have averted 24.8% to 
43.2%, if there was an increase in condom use. 
Discussion Increased condom use following the Avahan 
intervention is likely to have played a role in curbing the 
HIV epidemic in Mysore. In Belgaum, given the limitations 
in available data, this method cannot be used alone to 
decide if there has been an increase in condom use. 



INTRODUCTION 

Karnataka, in southern India, has one of the highest 
prevalences of HIV among Indian states, 1 and is 
among the states included in the 'Avahan' India HIV/ 
AIDS initiative established by the Bill & Melinda 
Gates Foundation in 2004. Working through state- 



level partners and local non-governmental organisa- 
tions (NGOs), the initiative is a multifaceted HIV 
preventive intervention promoting prevention 
strategies among high-risk groups, including distri- 
bution of condoms to and promotion of condom use 
among female sex workers (FSWs). 2-4 

Assessing the impact of the Avahan intervention 
is crucial, not only to ensure that Avahan achieves 
its goals of reducing HIV transmission in India, but 
also to inform future large-scale interventions. 5 
However, the dynamics of an HIV epidemic are 
complicated: declines in HIV prevalence can occur 
without any intervention effects, due to natural 
infection dynamics and, conversely, increases can 
occur even in the presence of an effective interven- 
tion. 6 7 In order to separate out the impact of the 
Avahan initiative from natural HIV transmission 
dynamics, and because of the difficulty and huge 
expense of implementing community randomised 
controlled trials, a tailor-made transmission 
dynamics model will be used as one important 
component of the impact assessment, the frame- 
work for which is described in Boily et al. B 

The final impact and cost effectiveness analysis, 
planned for 2011, will estimate the impact of the 
Avahan initiative on the HIV and sexually trans- 
mitted infection (STI) epidemics among high-risk 
groups and the broader population across more 
than 30 intervention sites. 5 8 This paper presents 
the results of an interim analysis focusing on FSWs 
and their clients in two districts where sufficient 
data are currently available, and which represent 
different HIV/STI epidemic trends and behavioural 
profiles: Belgaum Urban, a higher prevalence district 
in which Avahan was not the first but is now the 
only intervention; and Mysore Urban, a lower 
prevalence district in which Avahan was the first 
and remains the only intervention. 

The objective of this analysis is to determine 
whether the observed changes in HIV prevalence 
among FSWs could be solely due to natural disease 
dynamics rather than to a change in condom use 
following the intervention. We address this ques- 
tion by testing whether the null hypothesis (Ho) of 
stable (ie, no increase in) condom use following the 
initiation of the Avahan intervention in 2004 is less 
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likely than the alternative hypotheses (Hi) that condom use 
actually increased by as much as reported in FSW survey data, or 
(H 2 ) that there was a lower level of preintervention condom use 
and a sharper postintervention increase than indicated by FSW 
survey data, as suggested by records of condom availability. 9 

METHODS 

In order to model the HIV/STI epidemics in urban Belgaum and 
Mysore districts, an in-depth mathematical modelling analysis 
was undertaken with the key elements described below. 

Methodological framework of evaluation 

The structure of the evaluation follows the methodological 
framework designed and described in Boily et al. 8 In short, 
a mathematical model was used within a Bayesian framework, 
where available data were first used to specify what was known 
about each model parameter by defining a plausible range of 
values (the 'prior distribution'). This distribution was then 
randomly sampled repeatedly, to test the model with a large 
number of parameter combinations and identify those parameter 
sets (the 'posterior distribution') that agreed with the empirical 
HIV/STI prevalence data shown in table 1. 

One of the advantages of this approach is that the percentage 
of parameter sets producing a model fit to the prevalence data 
gives a measure of how compatible the prior parameter sets were 
with the observed prevalence data for the specified model. This 
approach was used here to test the three hypotheses described 
above concerning possible trends in condom use before and after 
the Avahan intervention. Finally the posterior distribution was 
used to simulate a control group without the intervention, 
which was then compared to the original scenario to produce 
estimates of intervention impact. 

Model structure and parameterisation 

A purpose-built deterministic compartmental model of HIV, 
herpes simplex virus (HSV)-2 and syphilis transmission was 
constructed and parameterised to incorporate key aspects and 
heterogeneities of the complex FSW and client structure specific 
to each setting. 10 In summary, the model consisted of an open 



Table 1 HIV/sexually transmitted infection (STI) prevalence data from 
integrated behavioural and biological assessment (IBBA) surveys in 
Belgaum and Mysore used to fit (FSW IBBA rounds 1 and 2, and client 
IBBA round 1) and validate (FSW IBBA rounds 2 and 3) the model 



Survey 



District and date Cl for P^^lence 
carried out 



HIV 



HSV-2 



Syphilis 



Fitting data 
FSW IBBA Mysore 
round 1 (August 2004) 
Belgaum 
(October 2005) 

Client IBBA 

Round 1 Mysore 

(October 2008) 
Belgaum 
(October 2007) 

Crossvalidation and fitting data 

FSW IBBA Mysore 

round 2 (December 2006) 
Belgaum 
(July 2008) 

Crossvalidation data 

FSW IBBA Mysore 

round 3 (April 2009) 



21 .9% to 30.3% 59.6% to 69.1% 21 .0% to 29.0% 
27.6% to 40.2% 78.6% to 89.1% 3.0% to 13.0% 



3.2% to 7.6% 8.0% to 33.0% 1 .3% to 4.6% 
3.6% to 8.8% 23.3% to 32.3% 2.0% to 6.5% 



19.1% to 29.5% 
22.2% to 32.5% 

8.11% to 14.1% 



population, growing at the rate described by census data, 1 1 stratified 
into 'high-risk' (FSWs and their clients) and 'low-risk' (the 
remaining population) groups. Based on work by Vickerman et at 
(see page 33), it was assumed in this model that there is 
no transmission between low-risk individuals. The high-risk 
groups (FSWs and their clients) were further stratified by variables 
showing strong associations with HIV prevalence and high-risk 
behaviour in analyses of data from serial cross-sectional surveys 
termed integrated behavioural and biological assessments 
(IBBAs). 10 12 

The model simulates the transmission of HIV/STIs between 
FSWs and their clients through commercial and longer-term 
non-commercial partnerships, and the bridging infections from 
clients to low-risk women through their non-commercial part- 
nerships, as shown in figure 1. The force of infection depends on 
the disease-specific infectivity, type and duration of partnership, 
frequency of sex acts for the type of partnership, condom use and 
STI coinfection status of both partners. The three parameters 
influencing HIV prevalence the most in a sensitivity analysis of 
a representative sample of the unrestricted model runs were the 
number of sex acts per partnership between FSWs and occasional 
clients, the male to female HIV transmission per sex act and the 
RR for increased transmissibility of primary HIV Further details 
of the model and equations are provided in the supplementary 
material. 

HIV was modelled with an initial short acute phase of high 
infectivity, followed by a long low-infectivity phase and a pre- 
AIDS phase of increased infectivity. 13 It was assumed that those 
with AIDS are chronically ill, and cease being sexually active. 
HSV-2 coinfection 14 and syphilis 15 were also modelled dynami- 
cally, with cofactors representing facilitation of HIV and HSV-2 
acquisition and transmission. 

The setting-specific behaviour-related, demographic-related 
and intervention-related prior model parameter distributions 
were derived from detailed serial cross-sectional behavioural and 
biological data (from IBBAs) collected from FSWs and their 
clients, and general population surveys (GPS) carried out as part 
of the Avahan monitoring programme, as well as from 
complementary data from other sources. The non-setting 
specific biological parameter ranges were based on reviews of 
relevant literature (tables 2 and 3). Prior distributions were 
assumed to be uniform across each parameter range; this was 
a 'conservative' choice to reflect the lack of precise knowledge 




FSW, female sex worker; HSV, herpes simplex virus. 



Figure 1 Outline of the main aspects of sexual behaviour structure and 
population movements between risk groups that are included in the 
model. HIV/sexually transmitted infection (STI) transmission through 
main (plain arrow) and casual (dotted arrow) partnerships is shown. 
Female sex workers and clients engage in commercial sex for a time 
before returning to the general population and being replaced from the 
general population (thick arrow). 



i34 



Sex Transm Infect 2010;86(Suppl 1):i33-i43. doi:1 0.1 136/sti.2009.038950 



Supplement 



and possibility of bias in the data used to estimate the 
parameters. 

The parameterisation of condom use before and after the 
intervention was data driven and specific to each hypothesis 
tested, resulting in three different prior condom use distribu- 
tions. The prior distribution of condom use for hypothesis Hi 
was based on estimates (with CIs) of site-specific condom use 
from 2001 onwards given by the analysis of Lowndes et al 75 That 
analysis used data from the IBBA round 2 surveys concerning 
time since starting to use condoms consistently with occasional 
clients to retrospectively 'reconstruct' the fraction of FSWs 
consistently using condoms each year from 2001 to the year of 
the survey. The null hypothesis H 0 assumed the same increase in 
condom use prior to 2004 as hypothesis Hi but no further 
increase after the start of the Avahan intervention in 2004. A 
further hypothesis (H 2 ) was based on HIV prevention 
programme data and other sources regarding the availability of 
condoms for FSWs 9 from 2004 to 2008 across Karnataka. This 
study suggested that the proportion of consistent condom users 
among FSWs in 2004 was much lower than suggested in 
hypothesis Hj due to low availability but that it rose steeply 
during 2006—2008 due to the Avahan intervention. The three 
different condom prior distributions are shown in figure 2, with 
the round 1 and 2 FSW IBBA point estimates shown for 
comparison. 

Parametric uncertainty analysis, model fitting and validation 

This step involved identifying combinations of model parame- 
ters that agreed with ('fit') the HIV/STI prevalence data from 
each site for each condom use hypothesis tested. Latin hypercube 
sampling (LHS) was used to sample several million parameter 
sets from all 'prior' parameter ranges for each hypothesis and 
district, with the sampling done so that the same set of param- 
eters, excluding the intervention parameters being tested, were 
used for the different hypotheses in each district. The model was 
run with each sampled parameter set to predict HIV/STI prev- 
alences. A parameter set was accepted as a fit if the model 
predictions simultaneously lay within the 95% CIs of the HIV, 
HSV-2 and syphilis prevalence data of the relevant FSW and 
client IBBAs, as described below. 

The HIV/STI prevalence data from FSWs and clients for each 
district were used in two stages: (1) as validation and (2) for 
hypothesis testing. In stage 1, the model was first fitted to just 
the round 1 FSW and client IBBA HIV, HSV-2 and syphilis 
prevalence data (see table 1), and the predictive value of these 
model fits was assessed against the HIV prevalence data not used 
directly at this stage of the fitting process, namely FSW HIV 
prevalence from the second round FSW IBBA in Mysore and 
Belgaum (table 1) and third round in Mysore (table 1), and FSW 
HIV prevalence by typology (street, brothel or home based) in 
rounds 1 and 2 (see supplementary material). In stage 2, the 
second round of IBBA HIV prevalence data were included in the 
fitting procedure along with the first round of FSW and client 
IBBA HIV/STI prevalence data (table 1) to produce a reduced set 
of model fits for hypothesis testing. 

Testing hypotheses of condom use 

The most likely condom use hypothesis was determined based 
on the relative proportions of fits to the first round FSW and 
client IBBA HIV/STI prevalence data and the second round FSW 
HIV prevalence data. For each hypothesis, the fraction of model 
fits from all parameter sets was compared to determine which 
hypothesis was more likely. This method, which is related to the 
version of Approximate Bayesian Computing as described by 



Weiss and von Haeseler, is also similar to the Bayesian methods 
used by Vickerman et al for determining the likelihood of 
different injecting drug user risk behaviour scenarios in 
Pakistan, 77 and Hallett et al for examining behaviour change in 
Zimbabwe. 78 

Impact analysis 

The posterior distribution of the most likely condom use 
hypothesis was used to simulate HIV prevalence and incidence 
over time following the start of the intervention, and in 
a matched control population. In the latter, each model fit was 
run again with condom use kept constant at 2004 levels (ie, at 
the start of the Avahan initiative). Comparing the impact and 
control scenarios, estimates of the impact of the increase in 
condom use following the start of the intervention among high- 
risk groups were derived for a 5-year period (2004—2009). These 
estimates derived from the posterior distribution include credi- 
bility intervals (CrI) generated by using all the parameter 
combinations that fit the prevalence data to reflect the uncer- 
tainty in parameter assumptions, since more than one set of 
parameters may produce an equally good fit to epidemiological 
trends. 

RESULTS 

Hypothesis testing 

The proportion of fits to each hypothesis is shown in table 4. In 
Mysore, the hypothesis H 2 , derived from analysing availability of 
condoms, gives the largest proportion of simulations which are 
consistent with prevalence data, and scenarios H 2 and Hi of 
increasing condom use following the introduction of Avahan are 
more likely than the null hypothesis H 0 of no condom increase 
after the start of Avahan, for which no model fits were obtained. 
In Belgaum, however, the null hypothesis H 0 of no increase use in 
condom use following Avahan and the alternative hypothesis Hi, 
based on the historical condom use reconstruction, cannot be 
distinguished by this method, although both are much more 
likely than the hypothesis H 2 . An alternative method to 
comparing the proportion of fits would be a likelihood-based 
approach, calculating the likelihood for each run; however this 
method gave comparable results (data not shown). 

Crossvalidation of results 

Table 4 shows that fits to only the round 1 FSW and client IBBA 
prevalence data (table 4) also generally fit the round 2 FSW IBBA 
HIV prevalence data well for all hypotheses (90% or more of 
posterior runs also fit FSW HIV round 2 prevalence). This then 
validates the model with respect to time trends in HIV (see also 
supplementary material). A third round of IBBA data was avail- 
able for Mysore only, and the model HIV prevalence projections 
were slightly higher (12.1% to 20.6% for H 2 ; 12.4% to 19.6% for 
Hi) than the third round IBBA data (8.1% to 14.1%). This is 
consistent with the fact that any further increases subsequent to 
the IBBA round 2 in condom use were not incorporated in the 
model. 

Lastly validation of the model by comparing its projections 
with FSW HIV prevalence by typology generally showed good 
agreement, although there was some residual heterogeneity, as 
might be expected (see supplementary material). 

Trends in HIV/STI prevalence and incidence 

The predicted prevalences over time for the model fits from the 
favoured hypotheses (H 0 and Hi for Belgaum and H 2 for Mysore) 
suggest that the HIV epidemics among FSWs and clients are 
declining in both districts (figure 3A— F). In Belgaum, it is likely 
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Table 2 Prior biological model input parameters sampled at the fitting stage to obtain the posterior parameter sets for Mysore and Belgaum for the 
Avahan impact model (all durations are in months) 



Types of model input Definition of model input 



Model inputs 



Reference for model input value 



Duration of 
infection stages 



Transmission 
probabilities 



Cofactors for HIV 



Average duration of Ng/Ct: 

Males 2-5 

Females 2-12 

Female sex workers (FSWs) 0.5—3 
Average duration of syphilis stages: 

Primary (no treatment) 1.51 

Secondary (no treatment) 3—4.5 

Primary and secondary stage (with treatment) 1—5 

Latent phase (including treatment) 2—24 

Time between potential recurrences 6 

Immune/resistant phase 12—60 
Average duration of HSV-2 stages: 

Primary stage 0.36—0.66 

Symptomatic recurrence 0.1—0.16 
Rate of HSV-2 symptomatic recurrences while: 

HIV negative 0.09-0.41 

HIV positive 

Average duration of HIV stages: 
Initial HIV high viraemia phase 
Between initial high viraemia and pre-AIDS 
Pre-AIDS high viraemia phase 
Probability of HIV transmission per sex act: 
Male to female: 
Female to male 

Sexual transmission multiplicative cofactor: 
Initial high viraemia phase 
Pre-AIDS high viraemia phase 
Probability of Ng/Ct transmission per sex act 
Probability of syphilis transmission per sex act (male to 
female): 

Ratio of transmission probabilities female to male:male to 
female 

Probability of HSV-2 transmission: 
Latent/asymptomatic shedding stage (male to female) 
RR male to female: female to male transmission 
Primary stage 

Symptomatic recurrence stage 



Reviewed in Korenrompef al. Also depends on level of STI 
treatment. 



Available data was reviewed by Boily ef al and 
Garnett ef a/ 8 15 



1-2XHIV negative rate 
4-6 

70-90.5 
6-18 

0.0006-0.0011 
0.0001-0.0014 

4.5-18.8 
4.5-11.9 
0.05-0.2 
0.1-0.3 

0.33-1.0 



From Cheong ef al, Corey ef al, Diamond ef al, Koelle ef al, 

Guinan ef al and Benedetti ef a/.' 7-24 

Using Corey ef al, Guinan ef al and Wald ef a/.' 8 20 23 25-27 



Corey ef al, Diamond ef al, Benedetti ef al, Lafferty ef al and 

Kim ef a/' 8 20 2 ' 24 28-30 

Using Schacker ef al and Conant ef a/ 31 32 

Based on Graver and Shivraj, and Kumarasamy ef a/ 33 34 



Reviewed in Holmes ef a/. 40 67 



Average Ng/Ct cofactor per sex act for increasing 
susceptibility to HIV 

Average syphilis cofactor per sex act for increasing 
susceptibility to HIV 

HSV-2 cofactor per sex act for increasing HIV infectivity: 
Primary and symptomatic recurrence phases 
Asymptomatic/latent phase 

HSV-2 cofactor per sex act for increasing HIV susceptibility 
Asymptomatic/latent phase 
Primary phase and symptomatic recurrence phase 



0.0005-0.002 

2-5 

2-6X6.7-25 times 
asymptomatic/latent transmission 
probability 
1-3X6.7-25 
times latent/asymptomatic 
transmission probability 
1.2-2.5 

2.1-3.3 



Reviewed by Holmes ef al, Hooper ef a/ 40 ' 
Reviewed by Garnett ef a/ 15 



Using Wald ef al and Corey ef al 



From Kim ef al, Wald ef al and Corey ef al 



Wald ef al, Kim ef al, and Mertz ef a/ 25 30 39-41 



1.27-2.57X1-2 
0.27-1.57X0.04-0.15X2-3 

1-4.75 

1.5—4.0 times asymptomatic 
cofactor 



Cofactors for HSV-2 HIV cofactor per partnership for increasing HSV-2 infectivity: 

Primary phase 1—2.5 

Asymptomatic/latent phase 2—4 

Symptomatic recurrence phase Same as primary 

Condom efficacies Condom efficacy per sex act for HIV 80% to 95% 

Condom efficacy per sex act for HSV-2 and syphilis 40% to 70% 

Condom efficacy per sex act for Ng/Ct 60% to 90% 



Using Rottingen ef a/ 13 



Using Schacker ef al, Nagot ef al, Zuckerman ef al, Baeten 
ef al, Celum ef al, Dunne ef al, Delany etal, Mbopi-Keou ef 
al, Augenbraunefa/, LeGoffefa/andSerwaddaefa/ 31 42-54 
and Quinn ef a/ 55 to convert from differences in HIV viral 
load 

Using Celum ef al, Corey ef al, Freeman ef al and Watson- 
Jones ef a/ 45 56-58 and converting to probability using 
Quinn ef a/ 55 



Pinkerton ef a/ 59 60 

Based on Wald etal, Oberle etal, Dobbins ef al, Obasi ef al 
and Huertaefa/ 37 6,-65 

Based on Holmes ef al, Hooper ef al, Austin ef al, Barlow, 
Joesof ef al, Sanchez ef al, Gaydos ef al, Niccolai ef al and 
Zenilman ef al 66-74 



Ct, Chlamydia trachomatis; HSV, herpes simplex virus; Neisseria gonorrhoeae; STI, sexually transmitted infection. 
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Table 3 Prior behavioural model 
Avahan impact model 



input parameters sampled at the fitting stage to obtain the posterior parameter sets for Mysore and Belgaum for the 



Types of model input 



Definition of model input 



Mysore 



Belgaum 



Demography 

Population size and demographic inputs 



Migration of female sex workers (FSWs) 



Sexual behaviour 

Long-term partnerships of FSWs and 
clients 



FSW sexual behaviour 



Client sexual behaviour 



Initial size of sexually active population 

Fraction of female general population 
sexually active 

Fraction of male general population 
sexually active 

Average time spent sexually active 

Entry rate into sexually active population 
per year 

Proportion of FSWs migrating 
Multiplicative cofactor increasing 
prevalence of HIV among clients in sites to 
which FSWs migrate 



278000 (M) 
268000 (F) 
75% to 95% 

Determined by available partnerships 

41.3 years (M) 

42.4 years (F) 
12000 (M) 
11000 (F) 
0.3-0.44 
1-2 



Percentage of clients currently married/cohabiting by duration: 
0-1 years and 2-4 years 68% 
5-9 years and 10+ years 68% 
Percentage of FSWs currently married/cohabiting by duration: 
0-1 years 38.2% to 57.2% 

2-4 years 39.9% to 54.9% 

5-9 years 28.6% to 50.4% 

10+ years 23.1% to 46/7% 

Percentage of FSWs currently married/cohabiting by typology: 



Home based 
Brothel based 
Street based 

Average frequency of sex acts with 
married/cohabiting partner for FSWs/ 
clients (per month) 
Duration of long-term partnerships 
(cohabiting/married) if FSW, years 
Duration of long-term partnerships 
between clients and low-risk females, 
years 

Average weekly frequency of clients for: 
Home-based FSWs duration 0—1 years 
Home-based FSWs duration 2—4 years 
Home-based FSWs duration 5—9 years 
Home-based FSWs duration 10+ years 
Brothel-based FSWs duration 0—1 years 
Brothel-based FSWs duration 2—4 years 
Brothel-based FSWs duration 5—9 years 
Brothel-based FSWs duration 10+ years 
Street-based FSWs duration 0-1 years 
Street-based FSWs duration 2—4 years 
Street-based FSWs duration 5—9 years 
Street-based FSWs duration 10+ years 
Number of sex acts with each client 
Average duration of sex work in months for: 
Home-based FSWs 
Brothel-based FSWs 
Street-based FSWs 
Number of FSWs visited/month if: 
Below median activity level 
Above median activity level 
Average duration of being client in months 
if below/above median activity level 



25.9% to 54.1% 
8.8% to 66.9% 
39.3% to 49.3% 
5.9-8.2 (clients) 
5.9-8.2 (FSWs) 

8.6-13.0 
20-30 



2.5- 5.1 
3.8-7.9 

2.8- 8.2 

4.6- 14.6 
15-25 
12-45 
4-15 
5.2-16.6 
6.5-9.2 
7.5-9.4 
7.2-9.8 
6.5-9.1 
1-3 

35-69 

49-189 

53-66 

0.8-1.2 

2.9- 3.8 

84-240/84-240 



268000 (M) 
257000(F) 
77% to 95% 

Determined by available partnerships 

41.5 years (M) 
42.2 years (F) 
12000 (M) 
11000 (F) 
0.05-0.20 
1-2 



43% to 59% 
83% to 93% 

9.8% to 40.3% 
11.6% to 39.9% 
21.5% to 48.6% 
5.3% to 20.7% 

1 1 .4% to 29.7% 
13.1% to 28.9% 
9.0% to 43.9% 

9.5- 11.7 (clients) 

5.6- 9.5 (FSWs) 

16.6-21.7 
20-30 



5.5-26.1 

1.8- 10.3 

5.7- 10.0 

7.0- 11.0 

6.3- 34.2 
13.5-26.2 
13.1-25.3 

8.9- 17.9 

4.8- 13.6 

1.1- 13.0 

2.4- 10.2 
3.8-10.2 
1-3 

160-224 

94-127 

138-251 

0.8-1.2 

2.32-2.78 

83-119/90-143 



Continued 
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Table 3 Continued 



Types of model input 



Definition of model input 



Mysore 



Belgaum 



Proportions of FSWs/clients by each 
stratification 



Condom use 

Condom use in main partnerships 
Condom use between FSWs and clients 



Dates 



Percentage of female population who are 
FSWs 

Percentage of male population who are 
clients 

Proportion of FSWs that are: 
Home based 
Brothel based 
Street based 

Proportion of male clients who visit FSWs: 
Below median level (ie, low activity 
clients) 

Above median level (ie, high activity 
clients) 

Average consistency of condom use 
between married/cohabiting partners per 
sex act. 

Fraction of sex acts with occasional clients for which a condom is used, by FSWs who: 
Report 'always' using 0.81—0.93 
Report often/sometimes using 0.54—0.67 
Report 'never' using 0.07—0.38 
Fraction of FSWs who are consistent condom users: 



0.3% to 1.4% 

Determined by number of FSW 
partnerships 

0.08-0.15 
0.00-0.01 
0.85-0.91 

0.5 
0.5 



4.3% to 10.3% 



0.2% to 1.1% 

Determined by number of FSW 
partnerships 

0.23-0.50 
0.31-0.61 
0.11-0.24 

0.58-0.69 
0.31-0.42 



5.2% to 12.7% 



At the start of the HIV epidemic in India 0-0.1 
Under H 0 and H,: 

At first time point 0.111-0.228 

At second time point 0.233-0.395 

Under H,, at time of IBBA R2 0.61 5-0.729 
Under H 2 : 

At first time point 0.1-0.220 

At second time point 0.321-0.416 

At time of IBBA R2 0.777-0.847 
Fraction of FSWs who report 'sometimes/ 0.245—0.367 
often' using condoms with occasional 
clients at IBBA R2 

Start of HIV epidemic in India 1976-1985 
Under H 0 and H,: 

At first time point 2002 

At second time point 2004 

At time of IBBA R2 2007 
Under H 2 : 

At first time point 2004 

At second time point 2006 

At time of IBBA R2 2008 



0.81-0.93 
0.54-0.67 
0.07-0.38 

0-0.1 

0.286-0.422 
0.759-0.855 
0.855-0.926 

0.1-0.220 
0.321-0.416 
0.777-0.847 
0.074-0.145 

1976-1985 

2002 
2005 
2008 

2004 
2006 
2008 



FSW, female sex worker; H, hypothesis; IBBA, integrated behavioural and biological assessments; R, round. 



that this decline began before Avahan started, around the year 
2000, while in Mysore the decline probably started shortly after 
Avahan, in 2004—2005. Modelled incidence and the prevalence 
for hypothesis Hi for Mysore show similar trends (see supple- 
mentary material). More details of the general characteristics of 
the model runs and fits are included in the supplementary 
material. 

In Mysore, the median relative decline in HIV prevalence in 
2009 compared to 2004 ('time decline') under H 2 was 36.2% 
(95% CI 21% to 51%) in FSWs (median 25.6%, 95% CI 7% to 
42% in clients), whereas the median relative decline in 2009 
prevalence compared to a simulated scenario without any 
increase in condom use since 2004 ('intervention decline') was 
44.7% (95% CI 35% to 52%) in FSWs (37.9%, 95% CI 30% to 45% 
in clients). In Belgaum, under H 1; the median time decline was 
47.6% (95% CI 39% to 58%) in FSWs (median 55.7%, 95% CI 
48% to 66% in clients) while the median intervention decline 
was actually lower at 23.1% (95% CI 17% to 30%) in FSWs 
(28.4%, 95% CI 22% to 36% in clients). Thus, the time decline in 



HIV prevalence slightly underestimates the intervention decline 
in Mysore for hypothesis H 2 and overestimates it in Belgaum for 
hypothesis Hi. 

Intervention impact among high-risk groups 

We estimated the impact of the increase in condom use 
following the start of Avahan using the most likely hypothesis 
as the simulated intervention group/population, and using the 
null hypothesis H 0 as the simulated matched control group/ 
population. In Mysore, the increase in condom use between 2004 
and 2009 may have prevented 31.3% to 47.4% of new HIV 
infections in FSWs and 32.7% to 47.2% in clients under 
hypothesis H 2 (figure 4 shows this in comparison to the impact 
projections for hypothesis Hi). Impact increases over time, with 
58.7% to 82.2% of new HIV infections prevented among FSWs in 
2008 and 64.1% to 85.1% among clients under hypothesis H 2 
(41.2% to 72.2% and 47.1% to 77.6% under Hi). 

In Belgaum, hypothesis H 0 , where condom use remains 
constant from 2004 onwards, suggests no impact. In 
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Figure 2 Prior parameter ranges for hypotheses H 0 — H 2 reflecting the 
proportion of female sex workers (FSWs) who are consistently using 
condoms for (a) Mysore and (b) Belgaum. The upper and lower lines 
show the range of the prior. Also shown is the data from the historic 
condom reconstruction, and mean and 95% CIs for the FSW integrated 
behavioural and biological assessment (IBBA) values. 

comparison, the Hi scenario projects that 24.8% to 43.2% of 
new FSW HIV infections and 30.9% to 53.1% of new client 
infections were prevented in Belgaum by the increase in condom 
use between 2004 and 2009 (figure 4). As for Mysore, impact 
increases over time in the Hi scenario, with 32.7% to 59.7% of 
new infections in FSWs and 42.0% to 72.6% in clients prevented 
in 2008. 



DISCUSSION 

With sufficient data to inform model parameters and validate 
model predictions, the use of mathematical models within 
a Bayesian framework enables the testing of hypotheses about 
the value or range of more uncertain parameters, 77 78 such as the 
evolution of condom use over time. In this analysis we have 
tested whether condom use is likely to have increased following 
the implementation of the Avahan intervention. In theory this 
could have been done more directly if unbiased baseline data on 
condom use had been collected; however, in practice it was not 
possible because such data could only be collected once the 
intervention was sufficiently well established to map and sample 
FSWs. 

The success rates obtained in the hypothesis-testing analysis 
support the conclusion that in Mysore condom use during 
commercial sex increased substantially after the beginning of the 
intervention. This increase most likely occurred shortly after the 
beginning of the intervention, as suggested by data on avail- 
ability of condoms. Although there was no significant decrease 
in the FSW IBBA prevalence data between rounds 1 and 2, the 
model suggests that, in the absence of an increase in condom use 
after 2004, prevalence would have increased, and so to obtain the 
required trend a significant increase in condom use was required. 
In Belgaum, the model results show that decreases in prevalence 
alone from the IBBA surveys are insufficient to distinguish 
between the hypotheses of increased condom use as reported by 
FSWs and constant condom use following the beginning of the 
intervention. Without information from other sources, it is not 
possible to definitely conclude that there has been an increase in 
condom use in commercial sex in Belgaum. The results also 
suggest that the level of condom use was higher in Belgaum 
than in Mysore prior to 2004. In Belgaum, a smaller-scale 
intervention was present for some years prior to 2004, 79 before it 
was scaled up by Avahan. Condom use by FSWs with clients may 
therefore have been higher at the start of the intervention than in 
Mysore, where there was no prior intervention. 

The model results indicate that HIV incidence and prevalence 
have declined over time. The decline occurred earlier in Belgaum 
than in Mysore because condom use became widespread earlier, 
but perhaps also because the epidemic matured earlier. In 
Mysore, the decline roughly coincides with the beginning of the 
intervention. This pattern of decrease in prevalence agrees with 
that seen by Vickerman et al 4 and by Boily et al 80 examining 
ANC data. While this study only considers high-risk individuals, 
the decline observed among FSWs and clients should translate to 
a slower decline in the general population in both districts, as 
shown in the analysis by Vickerman et al. 4 

The relative decline in HIV prevalence in Mysore in 2009 
compared to the simulated control group ('intervention decline') 
is higher than the decline over time since 2004 ('time decline'). 
This is because the epidemic was still growing in 2004. In 
Belgaum, under the hypothesis of increasing condom use after 



Table 4 Number of fits per million runs for each hypothesis (H) (fits to round (R)1 data only are shown for comparison to validate time trends of the 
model in HIV prevalence) 





Mysore 






Belgaum 








H 0 : condom use 
fixed since 2004 


H,: condom use as 
reconstructed from 
IBBA round 2 data 


H 2 : condom use follows 
condom availability trends 9 


H 0 : condom use 
fixed since 2004 


H,: condom use as 
reconstructed from 
IBBA round 2 data 


H 2 : condom use follows 
condom availability 
trends 9 


Fits to R1 FSW and 
client IBBAs 


0 


48.3 


86 


11.3 


9.7 


0 


Fits to R1+R2 FSW 
and client IBBAs 


0 


36 (S) 


82.7 (S) 


9.7 


8.3 |NS) 


0(S) 



FSW, female sex worker; IBBA, integrated behavioural and biological assessments; NS, not statistically different to H 0 at the 5% level using x 2 test; S, statistically significantly different to H 0 . 
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Figure 3 HIV prevalence over time for most likely hypothesis for (a) female sex workers (FSWs) and (b) clients in Mysore (H 2 ); (c) FSWs and 
(d) clients in Belgaum (Ho); and (e) FSWs and (f) clients in Belgaum (H,). Shown on the graphs are median, 25th and 75th percentiles, and the 95% 
credibility interval (lighter shaded area). For H, Belgaum and H 2 Mysore the prevalence of the simulated control groups is also shown (median, 25th and 
75th percentiles, and 95% credibility interval). Also shown is the integrated behavioural and biological assessment (IBBA) prevalence data. 



Avahan, the intervention decline is smaller than the time 
decline, as the epidemic there was already declining prior to 
Avahan. This highlights that time trends alone should not be 
used to evaluate an intervention, and that modelling results are 
important to complement prevalence data in carrying out an 
evaluation. 6 8 ^ 



Strengths and limitations 

The model described in this paper was specifically tailored to 
reflect the behavioural and biological heterogeneities present in 
these settings, with dynamical modelling of HSV-2 and syphilis 
in addition to HIV in each site. It was parameterised using data 
from multiple surveys among FSWs, clients and general 
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Figure 4 Fraction of new HIV infections averted over a 5-year period 
from 2004 in female sex workers (FSWs) and clients as a result of 
increases in condom use for (A) Mysore under hypotheses H, and H 2 and 
(B) Belgaum under hypotheses H,. Shown on the graph are the median 
(middle line in the box), 25th and 75th percentiles (box) and the 95% 
credibility interval (whiskers). Belgaum hypothesis H 0 corresponds to no 
impact, and is not shown. 



population, and captured key epidemiological trends in HIV, 
HSV-2 and syphilis by risk groups. This model also compared 
well to an independent model by Vickerman et al, 4 and it was 
also validated with respect to time trends in HIV and FSW HIV 
prevalences by typology. 

The framework described in Boily et al 8 was used to produce 
point estimates and credibility intervals of intervention impact 
by the creation of control groups that assumed no increase in 
condom use. This control group gives the impact due to the 
increase in condom use since the beginning of the intervention, 
and, in situations where condom use would have increased 
naturally even without an intervention to drive it, will over- 
estimate intervention impact. Conversely, in the absence of 
a likely increase in condom use following Avahan, as may have 
been the case in Belgaum, this choice of control group results in 
the model projecting no HIV infections averted, while it is 
possible that there would still have been a decline in condom use 
in this situation without an intervention to sustain it. Indeed, 
analysis of a slow, modest decrease in consistent condom use of 
10% over 5 years suggests that even such a small decline in 
condom use leads to a large impact (17.3% to 29.7% additional 
new FSW HIV infections between 2004—2009 due to the decrease 
in condom use), and thus the contribution of the intervention in 
sustaining condom use is important, although it is difficult to 
attribute directly. 

The precision of the estimates was limited by data availability. 
In both districts, the fitting success rate was modest because 
conservative (wide) prior ranges were used for most parameters 
to capture their uncertainty, and because the model had to fit 
seven separate prevalence outcomes. If more data were available, 
then it would be possible to use this additional information to 
reduce the width of these priors and so possibly identify a higher 
proportion of model fits. This would allow a stronger compar- 
ison between the condom use scenarios, and would better reflect 
uncertainty in the impact projections. HIV prevalence data for 
earlier in the epidemic would also enable us to better constrain 
the trajectory of the early epidemic and so rule out unlikely 
epidemics with high early prevalence. However, using 



a constraint such as FSW prevalence being <10% in 1985 did not 
change the impact estimates significantly (data not shown) . 

The picture is further complicated in Belgaum because the 
first IBBA survey was carried out 16 months after the start of 
the Avahan initiative, and so behavioural and epidemiological 
data reflected a situation in which Avahan was already present 
and where consistent condom use rates reported by FSW in the 
2005 IBBA were already very high, at 91%. In contrast, the 
survey in Mysore is likely to be more comparable to the pre- 
Avahan baseline as it took place 8 months after the start of 
Avahan. It is difficult to quantify the impact of a change in 
behaviour in situations where there was a pre-existing inter- 
vention and a long delay between the start of the intervention 
and the first IBBA survey, such as in Belgaum. Impact estimation 
is more uncertain in such settings because it is hard to define and 
simulate a control group reflecting how condom use would have 
decreased without the intervention, for which no data are 
available. An additional round of behavioural and HIV preva- 
lence data for FSWs in Belgaum could help to further inform 
impact estimates by giving extra information about how the 
epidemic there is changing. 

In Belgaum, the low number of fits from the model to the 
prevalence data may be largely explained by the difficulty of 
fitting the relatively high HIV prevalence in clients. This anal- 
ysis modelled the transmission of HIV/STIs among urban 
clients, since the IBBA was carried out in urban areas where the 
intervention started. However the size of the FSW population is 
higher in rural areas, and as city clients may also go to smaller 
towns and villages to visit FSWs, the urban client HIV preva- 
lence may be influenced by the large rural FSW population. 
Future data collection strategies for programme evaluation 
should be adapted to reflect the reality of how the epidemic is 
driven and the nature of the intervention. 

Conclusions 

Dynamical modelling, used within a Bayesian framework, is 
useful for testing hypotheses and providing informative and less 
subjective impact estimates on which to base decisions. This 
method of comparing proportions of fits of different scenarios is 
relatively straightforward to implement, as shown here and in 
another recent analysis, 77 and helps to further augment survey 
findings of behaviours that are prone to reporting biases, such as 
condom use. 

This analysis suggests that the HIV epidemics in Belgaum and 
Mysore are declining. In Belgaum, it is probable that the decline 
in HIV prevalence started before the beginning of the Avahan 
intervention, since condom use by FSWs with their clients was 
already high in this district due to pre-existing interventions. 
The model, using current prevalence data alone, is unable to 
distinguish between scenarios where condom use in commercial 
sex was sustained at existing levels and where it increased after 
Avahan. If there was a rise in condom use, as suggested by other 
data 9 and equally likely compared to the null hypothesis from 
this method, this increase would have averted 24.8% to 43.2% of 
HIV infections in FSWs and 30.9% to 53.1% of infections in 
clients between 2004 and 2009. If, however, Avahan only 
sustained the existing level of condom use attained by FSWs prior 
to Avahan, this is still important for controlling the epidemic, as 
even a small, gradual decline in condom use can greatly increase 
the number of new infections. This conclusion on impact, of 
including both scenarios, is conservative (which is deemed pref- 
erable), since the difference in proportion of fits between the 
hypotheses where condom use does or does not increase after 
2004 is small. However, combined with the additional data 
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Key messages 



► Transmission dynamics models used within a Bayesian 
framework are useful tools to help evaluate large-scale 
intervention in absence of a control group. Such a framework 
enables testing of different hypotheses and provides credibility 
intervals for estimates of infections averted. 

► Our modelling analysis suggests there has been a decrease in 
HIV prevalence in high-risk groups due to an increase in 
condom use since the beginning of the Avahan intervention in 
Mysore. 

► Due to data limitations, our analysis could not say conclusively 
whether the change in prevalence in Belgaum was due to an 
increase in condom use following the start of the intervention, 
or simply the result of the natural dynamics and increase in 




presented in Bradley ef al regarding the availability of condoms 
in Belgaum, it can be argued that it is likely that there has been an 
increase in condom use since the start of Avahan. 

In Mysore, the findings strongly suggest that condom use 
during commercial sex has increased since the introduction of 
the Avahan programme, and that this increase in condom use 
has averted 31.3% to 47.4% of HIV infections in FSWs and 
32.7% to 47.2% in clients in Mysore over the first 5 years. 
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